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Abstract 

This paper studies the Yang-Lee singularity of the 2-dimensional Ising 
model on the cylinder via transfer matrix and finite-size scaling tech- 
niques. These techniques enable a measurement of the 2-point and 3-point 
correlations and a comparison of a measurement of a corresponding uni- 
versal amplitude with a prediction for the amplitude from the (A4,Ai) 
minimal conformal field theory. 
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In 1985, Cardy [1] provided evidence that the Yang-Lee edge singularity [2] 
of the 2-dimensional (2D) Ising model could be described by the non-unitary 
(Ai,Ai) minimal conformal field theory (CFT) [3, 4] of the ADE classification 
[ ]. Cardy's identification provided several CFT predictions for the Yang-Lee 
edge singularity of the 2D Ising model. The predictions for the central charge, 
the exponent v, and the low energy excitation spectrum have been confirmed 
by numerical measurements on lattice spin models [ ] . 

Cardy also determined the forms of 2-point and 3-point correlations of the 
single primary field in the (A±,A\) minimal CFT [1]. These correlations define 
a universal amplitude, which is known as a structure constant [3, 7]. The pre- 
dictions of universal amplitudes are an important advance that CFT brought to 
the understanding of critical points. Since these predictions of CFT are based 
on nontrivial applications of crossing symmetries, they deserve to be experimen- 
tally tested [3]. 
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While lattice models have been used to numerically test some predictions 
of universal amplitudes in unitary minimal CFTs [3, 7, 8], no tests of such 
predictions have been done for the non-unitary CFTs. This article tests such 
a prediction for a non-unitary CFT by numerically measuring spin correlations 
at the Yang-Lee edge singularity of the 2D Ising model. The measurements 
enable a test of the CFT prediction of a universal amplitude in the non-unitary 
(Ai,A{) minimal CFT. 

As pointed out by Cardy [1], the non-unitary (A±,Ai) minimal CFT has a 
single primary field cf>(z,z), which has left and right conformal weights of -1/5 
and a scaling dimension x of -2/5. The 2-point and 3-point correlations of the 
primary field 4>{z, z) have the respective forms: 

G 4>4> (zi 1 z 1 ,z 2 ,z 2 ) = |(zi - z 2 \ 4/5 , (1) 

and 

G^(zi,zi,Z2,Z2,Z3,Z3,) = C\(zi - z 2 )(z 2 - z 3 )(z 3 - zi)| 2/5 - (2) 

With the Coulomb gas formalism [7], Cardy showed that the single structure 
constant. C, of the non-unitary (A^jAi) minimal CFT has the value [1]: 



/ rr(6/5)pr(i/5)r(2/5) 

V r(3/5)[r(4/5)]3 1 1 

In this article, numerical measurements at the Yang-Lee edge singularity check 
this CFT prediction. 

The numerical measurements were made for the 2D ferromagnetic Ising 
model whose Hamiltonian, H, is: 

M N 

H = - Y^-'^^S,.,.: + SijSi+u) + hS itj ]. (4) 



In the ferromagnetic Hamiltonian. the spin-spin coupling J is positive. For the 
Ising model, the Yang-Lee edge singularity occurs at temperatures above the 
critical temperature and for purely imaginary values of the magnetic field, h, 
i.e., h = iB with B real [2]. In particular, spin correlations were measured at a 
temperature, T, for which J/ksT = 0.1 

For the above Hamiltonian, the transfer matrix was used to evaluate the 
correlations. In particular, 2-spin and 3-spin correlations were evaluated on 
torii of length, M, and of various diameters, TV. In these evaluations, M was 
always much larger than N, i.e., M = 512 and N = 3 - 8. Thus, the measured 
correlations had the same distance behavior as correlations on an infinitely long 
cylinder when distances between spin fields were very small compared to M. 

From the numerical evaluations, finite-size scaling provided the tool for ex- 
tracting values of physical properties in the thermodynamic limit [9]. In partic- 
ular, the spin correlations are measured at special values of the purely imaginary 
magnetic field, h(N) = iByL{N). Each special value , Byl(N), satisfies the 
phenomenological renormalization group (PRG) equation for infinite cylinders 
of diameters (N — 1) and N: 

S(iB Y L(N),N-l) = Z(iB Y L(N),N) 

N - 1 N ' U 
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In the PRG equation, t;(iB, N) is the spin-spin correlation length on the infinite 
cylinder of diameter N when the magnetic field is iB. The PRG equation 
imposes that the spin-spin correlation length scales linearly with N as N — > oo. 
When evaluated at PRG values of the magnetic field, other physical quantities 
will scale with the width, N, to their values in the thermodynamic limit, i.e., 
near the Yang-Lee edge singularity [10, 6]. 

On an infinite cylinder of width N, CFT predicts that correlations will de- 
pend exponentially on the distances between the fields in the correlations when 
said distances are large compared the cylinder's diameter, N [11]. For a 2-point 
correlation, the correlation length is related to the scaling dimension, x, of the 
two conformal fields of the correlation. In particular, the 2-point correlation 
has the form exp(— 2irx(yi — y2)/N) when \y\ — 2/2 1 >> N. Here, y\ and j/2 are 
the positions of the two fields of the correlation along the axis of the infinite 
cylinder. For the 3-point correlation, the exponential behavior on the distances 
between the fields of the correlation is again fixed by the scaling dimensions, 
x's, of the various fields of the correlation. 

At the Yang-Lee edge singularity of the 2D Ising model, we used the am- 
plitudes of the 2-spin and 3-spin correlations, i.e., A ss and A sss , respectively, 
to evaluate the 3-spin structure constant. In particular, the measured values of 
the 3-spin structure constant, C(N), were obtained from the relation: 

_ A ass (iB YL (N)) 
° {N) ~ [A ss [iBy L {N)f' r { ) 

Here, A ss {iBYL{N)) and A sss (iBYL(N)), are the measured amplitudes of the 
2-spin and 3-spin correlations, respectively, at PRG values of the magnetic 
field. Similarly, we used the PRG measurements of the correlation length 
£{N) to measure the conformal dimension, x, of the spin field, i.e., x(N) = 
N/[2n£(iB Y L(N))]. The scaling behavior of these physical quantities with the 
cylinder's width, N, was used to obtain their values in thermodynamic limit, 
i.e., as N — > 00. 

Table 1 summarizes our transfer matrix results 1 for M = 512 and J/ksT = 

0.1. 



N 


Byl(N) 


x(N) 


\C(N)\ 


3 


0.184802 


0.353929 


1.80838 


4 


0.183348 


0.376870 


1.83711 


5 


0.183064 


0.385748 


1.85736 


6 


0.182982 


0.390108 


1.87054 


7 


0.182951 


0.392693 


1.87937 


8 


0.182946 


0.392911 


1.88633 


00 




0.398(2) 


1.923(13) 


CFT 




0.4 


1.9113 



Table 1: PRG Measurements of conformal dimension and structure constant. 



In these measurements of correlations, adjacent spin fields were separated 



1 While the B YL (N)'s for N = 3 - 7 were evaluated from the PRG equation, B Y l(N) for 
N = 8 was only approximately evaluated for numerical reasons. Byl{8) was determined from 
the Byl(NYs for N = 3 - 7 by assuming a leading finite-size scaling behavior. 
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by distances that were large compared to N and small compared to M. In 
particular, we measured correlations in which the distances between adjacent 
spin fields were between about AN and 0.06Af. For such distances, the CFT 
prediction of an exponential dependence on the distances between fields was 
obtained. In Table 1, the oo line shows values obtained by extrapolating our 
measurements to the thermodynamic limit, i.e., N — oo. The extrapolated 
values were obtained from nonlinear fits of the measured values of x(N) and 
|C(iV)| to functions of form f(N) = /(oo) + f\N~ a : i.e., to account for leading 
finite-size corrections. In each such extrapolation, the coefficient, /i, and the 
scaling power a were determined by finding the best fit to the measured data 
points. The best fits were found when the scaling corrections or x(N) and 
\C(N)\ had exponents, a, of about 2.4 and 1.2. In Table 1, the last line shows 
Cardy's CFT predictions for x and \C\ from the (A4, Ai) non-unitary minimal 
CFT as the model. 




Figure 1: The measured structure constant (squares) for cylinders of diameter 
N plotted against nonlinear fit for N = 4 — 8. 

Figure 1 plots our data for |C(iV)| and our best fit (line) based on a leading 
correction in N^ 1 ' 2 . In Figure 1, the black squares represent the measured 
C(7V)'s, the empty square represents the value of C(oo) obtained from the best 
fit to the measurements, and the black circle represents Cardy's CFT prediction 
for C(oo). 

The measured structure constants for cylinders of diameter 3-8 exhibits 
an approximate finite-size scaling correction of the form N~ 1,2 . The finite-size 
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scaling motivated fit matches the measured values of structure constant well 
except for the thinnest cylinder, N = 3. There, one is not surprised that a large 
N expansion does not work well. Finally, the dot shows the structure constant 
predicted by CFT. 

In conclusion, our finite-size scaling measurements for the 2D Ising model 
at the Yang-Lee edge singularity produce a value of the structure constant that 
agrees well with that of the (A4, A\) minimal CFT. These measurement confirm 
the CFT prediction of this universal amplitude at the Yang-Lee edge singularity 
of the 2D Ising model. 



References 

[1] J.L. Cardy, Phys. Rev. Lett. 54, 1354 (1985). 

[2] C.N. Yang and T.D. Lee, Phys. Rev. 87, 404 and 410 (1952). 

[3] A. A. Bclavin, A.M. Polyakov, and A.B. Zamolodchikov, Nucl. Phys. B241, 
333 (1984). 

[4] D. Friedan, Z. Qiu, and S. Shenker, Phys. Rev. Lett. 52, 1575 (1984), and 
Comm. Math. Phys. 107, 535 (1986). 

[5] A. Cappclli, C. Itzykson, and J.-B. Zuber, Nucl. Phys. B280, 445 (1987), 
and Comm. Math. Phys. 113, 1 (1987); A. Kato, Mod. Phys. Lett. A2, 
585 (1987); for a review see C. Itzykson and J.-M. Drouffe, Statistical Field 
Theory (Cambridge University Press, U.K., 1989) Ch. IX. 

[6] C. Itzykson, H. Salcur, and J.-B. Zuber, Europhys. Lett. 2, 91 (1986); K. 
Uzclac and R. Jullien, J. Phys. A14, L151 (1981); J.F. McCabc and T. 
Wydro, Int. J. of Mod. Phys. B20, 495 (2006). 

[7] Vl.S. Dotscnko and V.A. Fateev, Nucl. Phys. B240, 312 (1984), and B251, 
691 (1985). 

[8] G.T. Barkcma and J. McCabc, J. Stat. Phys. 84, 1067 (1996); J. McCabe 
and T. Wydro, Int. J. of Mod. Phys. A13, 1013 (1998); G. von Gehlcn, V. 
Rittenberg, and T. Vescan, J. Phys. A20, 2577 (1987). 

[9] M.E. Fisher and M.N. Barber, Phys. Rev. Lett. 28, 1516 (1972); for 
a review see M. Hcnkcl, Conformal Invariance and Critical Phenomena 
(Springer- Verlag, Germany, 1999) Ch. 3. 

[10] M.P. Nightingale, Physica A83, 561 (1976); and in Finite Size Scaling and 
Numerical Simulation of Statistical Systems, edited by Ed. V. Privman 
(World Scientific, Singapore, 1990) Ch. VII. 

[11] J.L. Cardy, J. Phys. A17, L385 (1984), and Nucl. Phys. B270, 186 (1986); 
for a review see M. Henkel, Conformal Invariance and Critical Phenomena 
(Springer- Verlag, Germany, 1999) Ch. 13. 



5 



